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ABSTRACT 

We present results of our Hubble Space Telescope Cycle 11 Survey for low-redshift (z < 1.65) damped 
Lya systems (DLAs) in the UV spectra of quasars selected from the Sloan Digital Sky Survey Early 
Data Release. These quasars have strong intervening Mgll-Fell systems which are known signatures 
of high column density neutral gas. In total, including our previous surveys, UV observations of Lya 
absorption in 197 Mgll systems with z < 1.65 and rest equivalent width (REW) W^ 2796 > 0.3 A have 
now been obtained. The main results are: (1) The success rate of identifying DLAs in a Mgll sample 
with W$ 2796 > 0.5 A and Fell W^ 2600 > 0.5 A is 36(±6)%, and increases to 42(±7)% for systems 
with W$ 2796 /W£ 2600 < 2 and Mgl W$ 2S52 > 0.1 A. (2) The mean HI column density of Mgll systems 
with 0.3 A < Wo A2796 < 0.6 A is (N(HI)) = (9.7 ± 2.2) x 10 18 cm~ 2 . For the larger REW systems 
in our sample, (N(HI)) = (3.5 ± 0.7) x 10 20 cm~ 2 . The mean HI column density remains constant 
with increasing REW for W^ 2796 > 0.6 A, but the fraction of Mgll systems containing DLAs is found 
to increase with increasing REW. (3) By combining our low-redshift results with results at higher 
rcdshift from Prochaska and Herbert-Fort and at z — from Ryan- Weber et al. we find that we can 
parameterize the DLA incidence per unit redshift as tidla(z) — no(l + z) 7 , where no = 0.044 ± 0.005 
and 7 = 1.27 ± 0.11. This parameterization is consistent with no evolution for z < 2 = 0.7, 
Qm = 0.3), but exhibits significant evolution for z > 2. (4) The cosmological neutral gas mass 
density due to DLAs is constant in the redshift interval 0.5 < z < 5.0 to within the uncertainties, 
Qdla ~ 1 x 10~ 3 . This is larger than Q. gas (z = 0) by a factor of w 2. (5) The slope of the HI 
column density distribution does not change significantly with redshift. However, the low redshift 
distribution is marginally flatter due to the higher fraction of high column density systems in our 
sample. (6) Finally, using the precision of Mgll survey statistics, we show that under the assumption 
of constant DLA fraction and HI column density suggested by our current sample, there may be 
evidence of a decreasing Qdla from z = 0.5 to z = 0. We discuss selection effects that might affect 
the results from our survey. We reiterate the conclusion of Hopkins, Rao, & Turnshek that very high 
columns of neutral gas might be missed by DLA surveys because of their very small cross sections, 
and therefore, that Qdla might not include the bulk of the neutral gas mass in the Universe. 
Subject headings: galaxies: evolution — galaxies: galaxy formation — quasars: absorption lines 



1. INTRODUCTION 

Recently, Fukugita & Peebles (2004) have summarized 
current measurements of the local mass-energy inven- 
tory. Of the local baryonic matter, about 6% is stars 
or their end states, about 4% is hot intracluster x-ray 
emitting gas, and somewhat less than 2% is neutral or 
molecular gas. The remainder of the baryonic matter is 
assumed to be in the form of a warm-hot intergalactic 
medium (WHIM), with properties similar to those dis- 
cussed by Cen & Ostriker (1999). However, importantly, 
processes in the neutral and molecular gas components 
most directly influence the formation of stars in galax- 
ies. Thus, the determination of empirical results on the 
distribution and cosmic evolution of neutral hydrogen 
gas is a key step in better understanding galaxy forma- 
tion. At present, there are two observational methods to 
study neutral hydrogen. Locally, the information is ob- 
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tained through radio observations of HI 21 cm emission. 
But at large distances (redshift z > 0.2) radio sensitiv- 
ity limitations require that the information be obtained 
through observations of Lya absorption in the spectra 
of background quasars. Intervening damped Lya (DLA) 
absorption-line systems in quasar spectra provide impor- 
tant non-local probes of the neutral gas content of the 
universe since they can, in principle, be tracked from the 
present epoch all the way back to the farthest detectable 
quasars. Since the first survey for DLAs nearly two 
decades ago (Wolfe, Turnshek, Smith, & Cohen 1986), 
it has been accepted that they contain the bulk of the 
neutral gas content of the universe. This first survey 
defined a DLA absorption-line system as an intervening 
gaseous HI region with neutral hydrogen column density 
N(HI) > 2xl0 20 cm~ 2 . The damping wings of the Voigt 
profile become prominent at column densities near 10 19 
cm~ 2 . Thus, even low-resolution spectra which are useful 
for the detection of Lya absorption lines with rest equiv- 
alent widths (REWs) > 10 A can be adopted to perform 
DLA searches, and subsequent studies have used this 
threshold to describe the statistics of DLAs (Lanzetta et 
al. 1991; Rao & Briggs 1993; Lanzetta, Wolfe, & Turn- 
shek 1995; Rao, Turnshek, & Briggs 1995, henceforth 
RTB95; Rao & Turnshek 2000, henceforth RT00; Storrie- 
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Lombardi & Wolfe 2000; Peroux et al. 2003; Prochaska 
& Herbert-Fort 2004). The N(HI) > 10 20 cm- 2 limit 
is believed to be the threshold above which the gas be- 
comes predominantly neutral and conducive for future 
star formation. 

The conclusion that DLA surveys identify the bulk of 
the neutral gas in the universe is based on three results or 
assumptions. First, integration of the DLA HI column 
density distribution shows that a relatively small frac- 
tion of the neutral gas is contributed by Lyman limit and 
sub-DLA absorption systems with 3 x 10 17 < N(HI) < 
2 x 10 20 cm" 2 , at least for z < 3.5 (Peroux et al. 2003, 
2005), and perhaps at all redshifts (Prochaska & Herbert- 
Fort 2004). Second, dust obscuration docs not cause 
DLA surveys to miss a large fraction of the neutral gas 
(Ellison et al. 2001, 2004). Third, the biases introduced 
by gas cross section selection are small. However, with 
regard to this last point, it is important to emphasize 
that the interception (or discovery) probability is the 
product of gas cross section times comoving absorber 
number density, and no DLAs with N(HI) > 8 x 10 21 
cm~ 2 have been discovered. Thus, the third assump- 
tion requires that rare systems with relatively low gas 
cross section and very high HI column density are ei- 
ther absent or have not been missed to the extent that 
the neutral gas mass density will be significantly under- 
estimated by quasar absorption line surveys. But this 
assumption might have to be reevaluated in order to ex- 
plain the discrepancy between the star formation history 
(SFH) of DLAs as inferred from their HI column den- 
sities and that determined from galaxies that trace the 
optical luminosity function (Hopkins, Rao, & Turnshck 
2005). We will, therefore, address the validity of this 
assumption later. 

Our main purpose in this paper is to present the re- 
sults of the most extensive survey for low-redshift DLAs 
to date. Since the Lya line falls in the UV for redshifts 
z < 1.65, Hubble Space Telescope (HST) spectroscopy 
is needed to detect and measure DLAs in this redshift 
regime that corresponds to the last w 70% of the age 
of the universe. Coupled with the fact that DLAs are 
rare, the scarcity of available HST time has meant that 
a good statistical description of the neutral gas content 
at low redshift is lacking. Now, with the Space Telescope 
Imaging Spectrograph (STIS) out of commission and the 
installation of the Cosmic Origins Spectrograph (COS) 
on HST only a remote possibility, further progress with 
UV spectroscopy seems unlikely, at least for the foresee- 
able future. 

To implement a low-redshift DLA survey with HST 
we have used an approach which differs from the conven- 
tional blind searches for quasar absorption lines. RTB95 
developed a method to determine the statistical proper- 
ties of low-redshift DLAs by bootstrapping from known 
Mgll absorption- line statistics. A similar approach was 
originally used by Briggs & Wolfe (1983) in an attempt 
to find 21 cm absorbers towards radio-loud quasars. It 
has been appreciated for some time that strong Mgll- 
Fell systems generally have HI column densities in excess 
of 10 19 cm -2 (e.g., Bergeron & Stasihska 1986). Thus, 
since all high-rcdshift DLAs are known to be accompa- 
nied by low-ionization metal-line absorption (e.g., Turn- 
shek et al. 1989, Lu et al. 1993, Wolfe et al. 1993, 
Lu & Wolfe 1994, Prochaska et al. 2003a and references 



therein), a UV spectroscopic survey for DLAs can be ac- 
complished efficiently if the search is restricted to quasars 
whose spectra have intervening low-ionization metal-line 
absorption. Since the MgIIAA2796, 2803 absorption dou- 
blet can be studied optically for redshifts z > 0.11, Mgll 
turns out to be an ideal tracer for low-redshift DLAs. 
If the incidence of metal lines is known, then the frac- 
tion of DLAs in the metal-line sample gives the incidence 
of DLAs. We further developed this method in RT00, 
and accomplished a three-fold increase in the number of 
low-redshift DLAs. We can now confidently use metal 
absorption-line properties as a predictor for the presence 
of DLAs. 

In this paper we present results from a sample of nearly 
200 Mgll systems with UV spectroscopy. Most of these 
data were obtained by us during the course of HST 
Guest Observer programs to make low-redshift surveys 
for DLAs. In principle, once DLAs are identified, follow- 
up observations can reveal details of a DLA's element 
abundances, kinematic environment, associated galaxy 
(i.e., a so-called DLA galaxy), star formation rate, tem- 
perature, density, ionization state, and size. For exam- 
ple, there is now clear evidence that the neutral gas phase 
element abundances are increasing from high to low red- 
shift (e.g., Prochaska et al. 2003b, Rao et al. 2005). 
There is a clear trend which indicates that DLAs resid- 
ing in regions exhibiting larger kinematic spread have 
higher element abundances (Nestor et al. 2003; Turn- 
shek et al. 2005). At low redshift (z < 1), it is now usu- 
ally possible to identify the DLA galaxy through imaging 
(e.g., Rao et al. 2003). At high redshift, high spec- 
tral resolution observations can be used to test dynami- 
cal models for DLA galaxies (Prochaska & Wolfe 1998). 
When the background quasar is radio loud, observations 
of 21 cm absorption provide important results on gas 
temperature (e.g., Kanckar & Chengalur 2003). Obser- 
vational constraints on physical conditions (temperature, 
density, ionization) also come from high-resolution spec- 
troscopy, and this has led to estimates of star formation 
rates in individual objects (Wolfe, Prochaska, & Gawiser 
2003). Estimates on the contribution of DLAs to the cos- 
mic SFH have also been made (Hopkins, Rao, & Turn- 
shek 2005). Finally, observations along multiple closely- 
spaced sightlines have led to estimates of DLA region 
sizes (Monier, Turnshek, & Rao 2005). Through such 
follow-up work, our knowledge of the characteristic prop- 
erties of the neutral gas component is steadily improving. 

Thus, our results serve two purposes. First, they pro- 
vide a comprehensive up-to-date list of more than 40 low- 
redshift (z < 1.65) DLAs suitable for follow-up studies. 
Second, they provide information on the distribution and 
cosmic evolution of neutral gas corresponding to the last 
?s 70% of the age of the universe. We discuss the Mgll 
sample in §2. The DLA sample is presented in §3, fol- 
lowed by statistical results derived from these systems in 
§4. Notably our study finds no evidence for evolution of 
the neutral gas mass of the universe between 0.5 < z < 5; 
at z = the neutral gas mass is now estimated to be a 
factor of « 2 lower. Moreover, at z < 2 there is no evi- 
dence for evolution in the product of absorber comoving 
number density and gas cross section, but at z > 2 their 
is clear evidence for an increase in this quantity in com- 
parison to no evolution models. A discussion of these 
and other new results is presented in §5. Conclusions are 
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summarized in §6. 

2. THE MGII SAMPLE 

The sample of Mgll lines used in our earlier DLA sur- 
veys (RTB95; RTOO) was culled from the literature. We 
observed 36 quasars that have 60 intervening Mgll ab- 
sorption systems with Mgll A2796 REWs W A2796 > 0.3 
A using HST-FOS in Cycle 6 (PID 6577). Twenty one of 
these Mgll systems fell in spectral regions with no flux 
because of intervening Lyman limit systems. Of the re- 
maining 39 systems, 9 were DLAs. With the addition of 
UV archival data, the total sample of Mgll systems with 
UV Lya information included 82 5 systems of which 12 
were DLAs. We found that all DLAs in this survey, with 
the exception of one, had Mgll W^ 2796 and Fell W$ 2600 
greater than 0.5 A. Based on this result, we conducted 
a similar survey of 54 Mgll systems in 37 quasars with 
HST-STIS in Cycle 9 (PID 8569). Most of these satisfied 
the strong Mgll- Fell criterion for DLAs. Twenty seven 
had useful UV spectra and four of these were DLAs. The 
DLA towards Q1629+120 was discovered in this survey 
and was reported in Rao et al. (2003). Results on the 
other systems from Cycle 9 are included in this paper. 

Further progress could only be made if the sample size 
was increased several fold. The Sloan Digital Sky Survey 
(SDSS) sample of quasars, which numbered in the thou- 
sands when this phase of our Mgll-DLA project began, 
presented an unprecedented leap in the number of avail- 
able survey quasars. The previous largest Mgll survey 
by Steidel & Sargent (1992; SS92) used a sample of 103 
quasars. The SDSS Early Data Release included nearly 
4000 quasars. Nestor (2004) used SDSS-EDR quasar 
spectra to search for Mgll systems with the aim of quan- 
tifying the statistical properties of a large Mgll sample 
(Nestor, Turnshek, & Rao 2005, henceforth NTR05) and 
to conduct follow-up work to search for DLAs. In Cycle 
11 (PID 9382), we targeted a sample of 83 Mgll systems 
with IU A2796 > 1 A in 75 SDSS quasars with SDSS mag- 
nitude g < 19. There were an additional 16 weaker Mgll 
systems observable in the same set of spectra. Overall, 
useful UV information was obtained for 88 systems, 25 
of which are DLAs. Given the large sample from which 
quasars could be selected for observation, we were able 
to minimize the occurrence of intervening Lyman limit 
absorption by restricting z em — z a b s to be small. Nev- 
ertheless, we were unable to observe the Lya line for 
11 Mgll systems either due to Lyman limits, an intrin- 
sic broad absorption line trough at the position of Lya 
absorption in one case, or due to the demise of STIS. 
Table 1 gives the details of the Mgll systems that have 
UV spectroscopic information. Our entire sample of 197 
Mgll systems is included. Details of the quasar are given 
in columns 1,2, and 3. Columns 4-8 give the absorption 

5 Of the 87 systems reported in Table 4 of RTOO, four sys- 
tems have been eliminated for reasons noted below, and one 
was reobserved in HST-Cycle 9. The z a i, a = 0.1602 system to- 
wards 0151+045 is a biased system because the galaxy-quasar pair 
was known prior to the identification of the Mgll system. The 
z abs = 0.213 system towards 1148+386 and the 0.1634 system to- 
wards 1704+608 were flagged as doubtful systems by Boisse et 
al. (1992). Also, on closer inspection, the IUE archival spec- 
trum of 1331+170 was inconclusive with regard to the Lya line 
of the 2 a f, s = 1.3284 system. Therefore, these four were eliminated 
from our current Mgll sample. The z a i s = 1.1725 system towards 
1421+330 is the one that was reobserved in Cycle 9. 



line information obtained either from the literature (col- 
umn 9 gives the reference) or from SDSS quasar spectral 
analysis (NTR05). Column 10 is the N(HI) measure- 
ment from the UV spectrum, column 11 is the selection 
criterion flag used to determine DLA statistics, and col- 
umn 12 gives the source of the UV spectrum. 

We now explain the selection criteria used to include 
Mgll systems in our sample. The total sample is di- 
vided into 4 sub-samples which essentially arose from 
the process of redefining and improving upon our se- 
lection process. Our first surveys, described in RTB95 
and RTOO, included strong Mgll systems from the lit- 
erature; the threshold REWs were chosen to match the 
Mgll sample of SS92 so that their statistical results could 
be used to determine the incidence of DLA systems. As 
demonstrated in RTOO, if the incidence of Mgll systems 
is known, then the fraction of DLAs in a Mgll sample 
gives the incidence of DLAs. We found that half of the 
Mgll systems in our sample with Wq 2796 > 0.5 A and 
Wq 2600 > 0.5 A were DLAs, and thus, modified our se- 
lection criteria to include a threshold Fell A2600 REW 
criterion as well. However, we retained the Mgll REW 
thresholds at 0.3 A, 0.6 A, and 1.0 A. 

The sub-samples are defined by the following criteria: 

1. W A2796 > 0.3 A; 

2. W A2796 > 0.6 A; 

3. Wo A2796 > 0.6 A and IU A2600 > 0.5 A; and 

4. IU A2796 > 1.0A and IU A2600 > 0.5 A. 

Sub-sample 1 includes all systems surveyed in RTOO, 
as well as additional systems that happened to fall along 
quasar sightlines that were targeted due to the presence 
of another stronger system from sub-samples 2, 3, or 4. 
Sub-samples 2 and 3 were mainly targeted for observa- 
tion in HST-Cycle 9, and sub-sample 4 includes systems 
found in SDSS-EDR spectra and observed in HST-Cycle 
11. A few systems from the SDSS-EDR sample have 
strong Mgll and Fell, but have Wo A2796 < l.oA; these 
belong in sub-sample 3. As we will see in §4.2, this 
classification is necessary for determining the incidence, 
it>dla{z), of the DLAs. 

3. THE DLAS 

Here we present Voigt profile fits to the new DLAs. 
These 28 systems, 3 of which were observed in HST- 
Cycle 9, are shown in Figure 1. The resulting column 
densities and errors are listed in Table 1. As is usually 
the case for high column density lines, the Lya forest 
populates DLA troughs making it inappropriate to use 
an automated routine such as least squares minimiza- 
tion to fit a Voigt profile to the data. Therefore, the best 
fit was estimated using the following procedure. Since 
the continuum fit is the largest source of uncertainty in 
determining N(HI), errors were determined by moving 
the continuum level by la above and below the best-fit 
continuum, rcnormalizing the spectrum, and refitting a 
Voigt profile (see RTOO). The differences between these 
values and N(HI) determined from the best-fit contin- 
uum are listed as the positive and negative errors in col- 
umn 10 of Table 1. 
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Fig. 1. — cont. 
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Fig. 1. — Voigt profile fits to the DLA lines. The quasar, Mgll 
z i,si an <J N(HI) are given in each panel. The dashed line is the 
best-fit Voigt profile and the dotted line is the la error array. 
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Fig. 2.— Plot of W A2796 vs. log N(HI). Filled circles are DLAs 
with N(HI) > 2 X 10 20 atoms cm . Arrows are upper limits in 
N(HI). Typical uncertainties are given by the error bars in the 
top left corner. 
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Fig. 3. — Distribution of MgIIA2796 rest equivalent widths, 
W/A2796. The shaded histogram represents systems that are DLAs. 
Note that there are no DLAs in the first bin, i.e., for Mgll 
W X2796 < 6 A. The fraction of DLAs increases with increasing 
W \2796 for W "A2796> Q g A. 
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4. STATISTICAL RESULTS 

4.1. Parameter Distributions and Correlations 

Since the systems in our sample were selected based on 
the rest equivalent width of Mgll A2796, measurements 
of Wo A2796 and N(HI) exist for all 197 systems. Mgll 
A2803, the weaker member of the doublet, was also mea- 
sured for all systems; measurements of the Fell A2600 
and Mgl A2852 lines were possible only for a subset. In 
this section, we explore correlations among metal line 
REWs and HI column density. Figure 2 is a plot of 
W A2796 versus log N (HI). We note that the upper left 
region of the figure is not populated, implying that sys- 
tems with Wq 2796 > 2.0 A always have HI column densi- 
ties N(HI) > 1 x 10 19 cm -2 . Figure 3 gives the distri- 
bution of Mgll W 7 ^ 2796 ; the DLAs form the shaded his- 
togram. It is noteworthy that there are no DLAs with 
Wq 2796 < 0.6 A. 6 In addition, the fraction of systems 
that are DLAs increases with increasing Wq 2796 . This is 
shown as a histogram in Figures 4 and 5; the y-axis on the 
left gives the fraction of DLAs as a function of Wq 2796 . 
We also plot the mean HI column density in each bin as 
solid circles with the scale shown on the right. Upper 
limits are assumed to be detections. 7 Figure 4 includes 
all observed Mgll systems and Figure 5 includes only the 
DLAs. The vertical error bars are standard deviations in 
the mean and are due to the spread of N(HI) values in 
each bin, and the horizontal error bars indicate bin size. 
For the Mgll systems there is a dramatic increase of a 



6 Only one known DLA has a lower metal-line REW. The 21 cm 
absorber at z = 0.692 towards 3C 286 has W$ 2796 = 0.39 A and 

1994). However, none of the 21 
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cm absorbers are included in our analysis because they are biased 
systems with respect to the determination of DLA statistics (see 
RT00). 

7 The two systems with high b values (the z a i, s = 1.6101 sys- 
tem towards 1329+412 and the z a i, a = 1.2528 system towards 
1821+107) are included in the histograms because it is clear that 
they are not DLAs. However, since their HI column density is not 
known, they are not included in the calculation of the mean column 
density. 




(A) 



Fig. 4. 



The histogram shows the fraction of Mgll systems that 
are DLAs as a function of Mgll Wq 2 " 79 ®, with the scale shown on 
the left axis. The solid circles are the logarithm of the mean HI 
column density in each bin. The scale is shown on the right. 



factor of ps 36 in the mean HI column density from the 
first to the second bin, beyond which (N(HI)) remains 
constant within the errors. In particular, Figure 4 shows 
that for our sample the probability of a Mgll system 
being a DLA is P « for W^ 2796 < 0.6 A and, assum- 
ing a linear dependence, P « 0.16 + 0.18(H / A2796 -0.6) 
for 0.6 < W n A2796 < 3.3 A. For systems with 0.3 A < 



W A2796 < 0.6 A, (N{HI)) = (9.7 ± 2.2) x 10 18 cm' 
and (N(HI)) = (3.5 ± 0.7) x 10 20 cm" 2 for systems with 
W A2796 > 0.6 A. Fi gure 5 shows a trend for decreasing 
DLA column density with Wq 2796 . The reasons for this 
are not obvious, but are likely to be due to small num- 
ber statistics (see Figure 2), a real physical effect, or a 
selection effect that is not yet understood (see Turnshek 
et al. 2005, §5.1). 
Since some of the higher Wq 2796 systems in the sam- 
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Fig. 5. — Same as Figure 4, but points are for DLAs only. 



pie have an FeIIA2600 selection criterion folded in, we 
repeat the above analysis here for the non-Fell selected 
part of the sample. Figure 6 is a W^ 2796 distribution for 
systems that do not include the Fell selection criterion. 
These are systems that belong to sub-samples 1 and 2. 
The DLAs form the shaded histogram. The first bin con- 
tains the same systems as in Figure 3. Even with this 
smaller sample, 111 systems compared to 197, we find 
that the fraction of DLAs increases with Wq 2796 . The 
mean HI column density for this sample is shown in Fig- 
ure 7. Since the number of systems in the higher Wq 2796 
bins is small, we bin the data differently from Figure 4 
and, for comparison, show the Figure 4 sample rebinned 
as well. We find that the Fell selection has no effect on 
the mean column density as a function of Wq 2796 . Con- 
sistent with our larger sample, for systems with Wq 2796 > 
0.6 A we find (N(HI)) = (3.40 ± 1.25) x 10 20 cm- 2 . One 
might expect that since the fraction of DLAs increases 
with increasing Wq 2796 and that the Fell selection pri- 
marily affects higher Wq 2796 systems, the mean HI col- 
umn density should be higher in the Fell selected sample. 
However, again, this may be offset by the fact that the 
mean DLA HI column density decreases with increasing 
Wq 2796 , thus keeping the mean HI column density of Fell 
and non-Fell selected samples indistinguishable. 

Figure 8 is a plot of W^ 2796 vs. W^ 2600 for systems 
with measured values of Wq 2600 , including upper limits. 
In RT00 we found that 50(±16)% of the 20 systems (ex- 
cluding upper limits and 21 cm absorbers) with Wq 2796 
> 0.5 A and Wo A260 ° > 0.5 A are DLAs. Now, with 
the expanded sample that includes 106 systems in this 
regime, we find that 36(±6)% are DLAs. The dashed line 
is a least-squares fit with slope b = 1.36 ± 0.08 and inter- 
cept a = 0.24 ± 0.06. It was determined using the BCES 
estimator of Akritas & Bershady (1996), assuming intrin- 
sic scatter but uncorrelated errors in Wq 2796 and Wq 2600 . 
Upper limits were not used for the fit. We note that 
DLAs do not populate the top left region of the diagram 
where the Wq 279Q to Wq 2600 ratio is > 2. In fact, if the 
sample is restricted to systems with Wq 2796 /Wq 2600 < 2, 
all but one of the DLAs in Figure 8 are retained, the out- 
liers in the top left region are excluded as are most sys- 
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Fig. 6. — Same as Figure 3 but for non-Fell selected systems. 
These form sub-samples 1 and 2. The shaded histogram shows the 
DLAs. As in Figure 3, the fraction of DLAs in bins with Mgll 
W/A2796> 0.6 A increases with increasing Wq 27ss . 



S5 

V 





i i i i i i i i i i 









w„ 



(A) 



Fig. 7. — Logarithm of the mean HI column density of absorbers 
as a function of Mgll Wq 27 ® 6 . The red solid squares are for non- 
Fell selected systems, i.e., for sub-samples 1 and 2. The black 
solid circles are the same data shown in Fig. 4, but rebinned to 
match the binning of the sub-sample without Fell selection. The 
data points in the first bin are identical but have been displaced 
for clarity. The Fell selection has no effect on the mean HI column 
density as a function of Mgll Wq 2796 . 



terns in the lower left corner of the plot. Figure 9 shows 
this truncated sample; the slope of the least-squares fit 
does not change significantly. We find b — 1.43 ± 0.08 
and a — 0.01 ± 0.08 for this definition of the sample. 
The only DLA that has been eliminated is the one with 
the smallest value of Wq 2600 . However, given the mea- 
surement errors for this system, its W / A2796 /W / ( ^ 2600 ratio 
is within la of 2. The implication is that a system with 
metal line ratio Wq 2796 /Wq 2600 > 2 has nearly zero prob- 
ability of being a DLA. For this truncated sample with 



W^ 27 "" 



IW> 



A2600 



< 2, but no restrictions on the individ- 



ual values of VF A2796 or W A2600 , 38(±6)% are DLAs. 
In addition, all known 21 cm absorbers, including the 
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(A) 



Fig. 8.— Plot of W A2796 vs. W A260 ° for all Mgll systems that 
have measured values of Wq 2600 . Filled circles are DLAs. Typical 
error bars are shown at lower right. The dashed line is the best fit 
linear correlation described in the text with slope b = 1.36. 
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Fig. 10.— Plot of W A2796 /W A260 ° vs. N(HI). Systems with 
Wq 2796 /Wq 2600 > 5 are not shown for clarity; all of these have 
log N (HI) < 19.6. The DLAs (filled circles) are confined to the 
region of the plot where 1 <IV A2796 /W A260 ° < 2. 
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Fig. 9.— Plot of W, 



W, 



A2796 



X2600 



(A) 



W A2600 



for systems with 

VF A2796 /W A2600 < 2. Filled circles arc DLAs. Non-DLAs in the 
upper left and lower left region of the diagram have been elimi- 
nated and the correlation is tighter than what is seen in Figure 
8. Also, 38% of all systems are DLAs regardless of the values of 
Wq 2796 and Wq 2600 ; we believe this to be a more robust predictor 
of the presence of DLAs in Mgll-Fell systems. The slope of the 
best fit linear correlation is b = 1.43. 



z = 0.692 system towards 3C 286 mentioned above, have 
W A279 7^o A2600 < 2. Thus, the VF A2796 /W X260a ratio 
provides a more robust predictor of the presence of a 
DLA. 

This result is shown more dramatically in Figures 
10 and 11. The ratio 1F A2796 /VF A260 ° is plotted as a 
function of N(HI) in Figure 10. Ratios above 5 are 
not shown for clarity. These are mainly confined to 
log N (HI) < 19.2 with only one system above this col- 
umn density at log N (HI) = 19.6. The DLAs populate 
the region of the plot where 1 <VF A2796 /VF A2600 < 2; the 
two outliers lie within lcr of this range. A plot of the 
ratio T^ A2796 /W A260 ° vs. VF A2852 for systems with mea- 




Fig. 11. 



0.5 



Plot of W A2796 



1.5 



(A) 

vs. Mgl 



Filled 



circles are DLAs. Typical error bars are shown in the lower 
right corner. Again, the DLAs are confined to the region where 
1 <W A2796 /W A2600 < 2, but span almost the entire range of 
WA2852 (W A2852 > 0.1 A). 



sured values of W 7 ^ 2852 , including upper limits, is shown 
in Figure 11. Again, the DLAs are confined to the region 
where 1 <W A2796 /W^ 2600 < 2, but span almost the entire 
range of W A2852 (W A2852 > 0.1 A). The two systems out- 
side the range 1 <W$ 2796 /W$ 2600 < 2 from Figure 10 do 
not have information on Wq 52 . Of the systems with 
measured values of TL^ 2852 , 32 of the 77 systems with 
VF A2796 /W A2600 < 2 and VF A2852 > 0.1, i.c, 42(±7)%, are 
DLAs. The other 9 DLAs cither do not have measured 
values of TL^ 2852 , or have high upper limits due to poor 
data quality. We also find that 9 out of the 11 systems 
with VF A2852 > 0.8A are DLAs. We note that systems 
with W A2796 /W A2600 > 2 are likely to have low values of 

W Q A2852_ 

For completeness, we also plot W^ 2600 vs. log N (HI) 
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Fig. 12.— Plot of W A2600 vs. log N (HI). Arrows indicate 
upper limits. Filled circles are DLAs. Typical uncertainties are 
given by the error bars in the top left corner. 
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Fig. 13. — Plot of W A2852 vs. log N(HI). Arrows indicate upper 
limits. Filled circles are DLAs. Typical uncertainties are given by 
the error bars in the top left corner. 



in Figure 12 and W A2852 vs. log N (HI) in Figure 13. 
There is no obvious trend in these distributions except 
for the fact that the upper left regions of the plots are 
not populated. There are no high REW, low HI column 
density systems. This is not a selection effect since col- 
umn densities as low as 10 18 cm" 2 can often be easily 
measured. This implies that systems with Wq 2W0 > 1 
A or Wq 2S52 > 0.5 A generally have HI column densities 
N(HI) > 10 19 cm' 2 . Below this fairly sharp boundary, 
metal-line REWs span all values of HI column density. 

4.1.1. Discussion 

How can these trends be interpreted? Apart from the 
upper envelopes in Figures 2, 12, and 13, there is no 
other simple correlation between metal-line rest equiva- 
lent width and HI column density. Since the metal lines 
are saturated, the rest equivalent width is more a mea- 
sure of velocity spread, not column density. High reso- 
lution observations of Mgll absorption lines have shown 



that the stronger systems break up into many compo- 
nents (e.g., Churchill, Vogt, & Charlton 2003), and span 
velocity intervals of up to 400 km/s. Turnshck et al. 
(2005) show line equivalent widths in velocity units of 
> 800 km/s in the strongest systems found in the SDSS. 
These highest equivalent width systems may arise in 
galaxy groups; however, the more common systems like 
those in our DLA survey are more likely to arise in clouds 
that are bound in galaxy-sized potentials. A DLA is ob- 
served if at least one of the clouds along the sightline 
happens to be cold (less than a 100 K), and with a ve- 
locity dispersion of a few 10s of km/s. A simple inter- 
pretation is that the greater the number of clouds along 
the sightline, the higher the probability of encountering 
a DLA. This would explain the higher fraction of DLAs 
among large Wq 2796 systems and the lack of a correla- 
tion between Wq 2796 and N(HI) other than the upper 
envelopes in Figures 2, 12, and 13. Only rarely would a 
sightline intersect a single cloud resulting in small W^ 2796 
and high N(HI), as in the 3C 286 system described in 
§4.1. This probabilistic approach to explain metal-line 
and HI strengths in high- N (HI) absorbers was also pro- 
posed by Briggs & Wolfe (1983) to explain their Mgll 
survey for 21 cm absorbers. They proposed a two-phase 
model where the 21 cm absorption is produced in galaxy 
disks, and the metal-line components that do not pro- 
duce 21 cm absorption are produced in galactic halos. 
However, this multi-component/cloud model is likely to 
be valid in any gas-rich galaxy, as is evidenced by DLA 
galaxy imaging studies (Le Brun et al. 1997; Rao & 
Turnshek 1998; Turnshek et al. 2001; Rao et al. 2003; 
Turnshek et al. 2004). The disk models of Prochaska 
& Wolfe (1997) and the Haehnelt, Steinmetz, & Rauch 
(1998) models of infalling and merging clouds could re- 
produce these observations equally well. In other words, 
DLAs arise in pockets of cold gas embedded within warm 
diffuse gas or gas clouds in any bound system. 

Twenty one cm observations of low-rcdshift DLAs also 
reveal some cloud structure. For example, the z = 0.313 
system towards PKS 1127—145 shows 5 components and 
the z = 0.394 system towards B0248+430 is resolved 
into 3 components (Lane 2000; Lane & Briggs 2001; 
Kanekar & Chengalur 2001). Since the Mgll line for 
these systems has not been observed at a resolution as 
high as the 21 cm observations, a one-to-one correspon- 
dence between the metal- line and 21 cm clouds cannot 
be drawn. In other instances, both warm and cold gas 
have been detected in a 21 cm DLA; Lane, Briggs, & 
Smette (2000) find that two-thirds of the column density 
in the z = 0.0912 DLA towards B0738+313 is contained 
in warm phase gas, and the rest is contained in two nar- 
row components. The z = 0.2212 absorber towards the 
same quasar was also found to exhibit similar charac- 
teristics (Kanekar, Ghosh, & Chengalur 2001). In each 
of these cases, the line of sight probably intersects two 
cold clouds in addition to warm diffuse gas spread over 
a wider range of velocities that can be detected only in 
21 cm observations of very high sensitivity. There are 
also several instances of DLAs not being detected at 21 
cm (Kanekar & Chengalur 2003). High spin tempera- 
tures (T s > 1000 K) corresponding to warm diffuse gas 
and/or covering factors less than unity towards extended 
quasar radio components have been suggested as possible 
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explanations (Kanekar & Chengalur 2003; Curran et al. 
2005). 

Clearly, a wide variety of cloud properties and their 
combinations are responsible for the observed properties 
of Mgll, DLA, and 21 cm absorption lines. Large simu- 
lations of galaxy sightlines with varying cloud properties 
that reproduce the metal-line versus DLA correlations 
shown in Figures 2-13 would be an important next step 
towards improving our understanding of these absorption 
line systems. The simulations should not only be able to 
reproduce the frequency of occurrence of DLAs in Mgll 
systems, but also the number density evolution of Mgll 
systems and DLAs. Moreover, further analysis on large 
data sets might enable us to predict the occurrence of 
DLAs among metal-line systems and determine their HI 
column densities to some degree of accuracy, but this is a 
project for future study. For the remainder of this paper 
we discuss the statistical properties of neutral gas in the 
low redshift universe as derived from the expanded HST 
sample. 

4.2. Redshift number density u^la 

The redshift number density of DLAs, udla, some- 
times written as dn/dz, can be determined using the 
equation 

tidla{z) = n(z)n Mg ii(z), (1) 

where n(z) is the fraction of DLAs in a Mgll sample 
as a function of redshift and nMgii(z) is the redshift 
number density of Mgll systems. Since our Mgll sample 
was assembled under various selection criteria (see §2), 
nMgii(z) needs to be evaluated carefully. We can express 
nMgii(z) for our sample as 

n Mg ii{z) = j^J2 WtnM 9 H ,( z )i ( 2 ) 

i 

where the sum is over all 197 systems, Wi is a weighting 
factor that depends on the i th system's selection criterion 
for being included in the survey, and nugii i (z) is the i th 
system's dn/dz value calculated using the parametriza- 
tion derived in the Appendix of NTR05: 

dn/dz = N* (1 + z) a e-^( 1+z ^\ (3) 

where we have retained the notation given in NTR05 
and W = W^ 2796 . N*, W* , a, and are constants. 
This expression is an integral over all Wq 2796 greater 
than Wo- For our calculation, Wo is different for each of 
four sub-samples that comprise our total sample (see §2 
and Table 1). Thus, for example, a system that belongs 
to sub-sample 1 has a REW threshold Wo = 0.3 A in 
Equation 3 and weight u>i = 1 in Equation 2, while a 
system in sub-sample 2 has Wo = 0.6 and to* = 1. This 
is because sub-samples 1 and 2 are purely Mgll-sclcctcd 
samples with no regard to the strength or presence of 
the Fell A2600 line. On the other hand, a system that 
belongs to sub-sample 3 has Wo = 0.6 in Equation 3 and 
weight Wi — 0.54 in Equation 2. This is because an Fell 
A2600 criterion was used to select the system in addi- 
tion to W* 2796 , and 54% of the 1,130 W A2796 > 0.6 A 
systems in the Mgll survey of NTR05 have Wo A2600 > 0.5 
A. Similarly, for systems in sub-sample 4, Wo = 1.0 and 
Wi = 0.72. ' In this case, 72% of the 781 W A2796 > 1.0 A 
systems in NTR05 have W A2600 > 0.5 A. For sub-samples 



3 and 4 we have assumed that the fraction of Mgll sys- 
tems that are also strong Fell systems is independent of 
redshift. 

The Mgll doublet moves out of the SDSS spectro- 
scopic range for rcdshifts z < 0.36. In order to ex- 
tend Mgll statistics to lower redshifts, we conducted a 
survey of quasars with the Multiple Mirror Telescope 
(MMT) on Mount Hopkins, AZ (Nestor, 2004). These 
results, which will be presented in a forthcoming paper 
(Nestor, Turnshek, & Rao 2006, in preparation), were 
used to determine DLA statistics for the redshift range 
0.11 < z < 0.36 using the same procedure described 
above. Of the 11 systems from our sample in this red- 
shift range, nine belong in sub-sample 1 and two are in 
sub-sample 2. 

Figure 14 shows the results for noLA(z) at low red- 
shift split into two redshift bins (solid squares). We find 
18 DLAs in 104 Mgll systems in the redshift interval 
0.11 < z < 0.9 with n DLA {z = 0.609) = 0.079±0.019 and 
23 DLAs in 94 Mgll systems in the redshift interval 0.9 < 
z < 1.65 with n DLA (z = 1.219) = 0.120 ± 0.025. We did 
not find it necessary to apply a Malmquist bias correc- 
tion for the number of systems with N(HI) > 2 x 10 20 
cm~ 2 (as was done in RT00) because the sample con- 
tains an equal number of systems within la above and 
below this threshold value. Standard error propagation 
procedures were used to determine uncertainties. The 
points arc plotted at the mean redshift of the Mgll sam- 
ples. The high-redshift data points are from Prochaska 

6 Herbert-Fort (2004) and the z = point was esti- 
mated by Zwaan et al. (2005a) from a WSRT survey 
of HI in the local universe. The solid curve is a no- 
evolution curve in the standard ACDM cosmology that 
we refer to as the "737" cosmology where (/i,Qm,^a) = 
(0.7, 0.3, 0.7). This curve, which has been normalized 
at the z = data point, shows that the comoving cross 
section for DLA absorption declined rapidly by a factor 
of w 2 until z « 2 and has remained constant since then. 
This behavior might be a consequence of what has been 
observed in other studies of galaxy evolution, namely, 
that today's galaxies were in place by z w 1 and are a 
consequence of rapid merger and/or collapse events that 
occurred prior to this epoch. 

It has been customary in quasar absorption line stud- 
ies to plot the logarithm of the redshift number density 
in order to illustrate its power law dependence with red- 
shift, i.e., n(z) = no(l + z) 1 . In A = cosmologies, the 
exponent is a measure of evolution. For example, 7 = 1 
for qo = or 7 = 0.5 for qo = 0.5 implies no intrin- 
sic evolution of the absorbers. Any significant departure 
from these values for 7 was considered as evidence for 
evolution in the product of the comoving number den- 
sity and cross section of absorbers. We plot log Udla{z) 
as a function of log(l + z) in Figure 15. The straight 
line is the power law fit to the data points with slope 

7 = 1.27 ± 0.11, and the curve is the same no-evolution 
function shown in Figure 14. Thus, in the past, the obser- 
vations would have been interpreted as being consistent 
with the DLA absorbers undergoing no intrinsic evolu- 
tion in a go = universe, and marginally consistent with 
evolution in a q = 0.5 universe. With the now widely 
accepted concordance cosmology, the interpretation has 
changed quite dramatically; as noted above, the nature 
of the evolution is redshift dependent. 
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Fig. 14. — Plot of noLA(z) versus redshift. The new low-redshift 
data points are shown as filled squares. The high-redshift points 
(filled circles) are from Prochaska & Herbert-Fort (2004) and the 
z = data point is from an analysis of local HI using the WSRT 
(Zwaan et al. 2005a). The solid line is a no-evolution curve in 
a standard "737" ACDM cosmology with (h,U M ,U A ) = (0.7, 0.3, 
0.7) normalized at the 2 = data point. This curve implies that 
the comoving cross section for absorption declined rapidly by a 
factor of Ri 2 until z ss 2 and has remained constant since then. 
This is consistent with the idea that today's structures have been 
in place since z 1 and are a consequence of merger events that 
occurred prior to this epoch. 



a 

60 




0.3 



0.4 



0.6 



log(l + z) 



Fig. 15. — Plot of \ogriDLA{z) as a function of log(l + z). The 
straight line is the power law fit to the data points with slope 
7 = 1.27 ± 0.11, and the curve is the no-evolution function shown 
in Figure 14. 



Further implications of this evolution are discussed 
in §5 along with inferences drawn from the evolution 
in Qdla (§4.3) and the HI column density distribution 
(§4-4). 

4.3. Cosmological mass density Qdla 

We can determine VLdla from the DLA column densi- 
ties listed in Table 1 and udla(z) via the expression 



^dla(z) 



fj,m H H 
cp c 



n DLA {z) (N(HI)) 



E(z) 
[l + zf 



(4) 



Fig. 16. — Cosmological mass density of neutral gas, Qdla, 
as a function of redshift. The filled squares are the new low- 
redshift data points. The high-redshift points (filled circles) are 
from Prochaska & Herbert-Fort (2004) and the open circle at z = 
is from Zwaan et al. (2005b). The open square at z = is the mass 
density in stars estimated by Panter, Heavens, & Jimenez (2004) 
from the SDSS. While the statistics have improved considerably, 
our basic conclusion from RT00 has remained unchanged, namely, 
that the cosmological mass density of neutral gas has remained 
constant from z sj 5 to z si 0.5. 



where 



E(z) 



H(z) 



p Af (i+z) 3 +(i-r! A/ -r! A )(i+z) 2 +fi A ] 1 / 2 . 

(5) 



Again, the "737" cosmology has been used in the calcu- 
lation of VLdla- Also, (J, = 1.3 corrects for a neutral gas 
composition of 75% H and 25% He by mass, mn is the 
mass of the hydrogen atom, p c is the critical mass den- 
sity of the universe, and (N(HI)) is the mean HI column 
density of DLAs in each bin. 

In contrast to the redshift number density evolution 
shown in Figure 14, we find that £Idla has remained 
constant from z — 5 to z = 0.5 to within the uncertain- 
ties. Figure 16 shows the new results as solid squares. 
Specifically, for the redshift range 0.11 < z < 0.90, 
we find (N(HI)) = (1.27 ± 0.36) x 10 21 cm" 2 and 
fl DLA (z = 0.609) = (9.7±3.6) x 10~ 4 , and for the range 
0.90 < z < 1.65 we get (N(HI)) = (1.07 ± 0.23) x 10 21 



cm" 2 and Qdla(z = 1.219) = (9.4 ± 2. 
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uncertainties have been reduced considerably in compar- 
ison to our results in RT00. The reasons for this are 
two fold. First, the uncertainty in riMgii has been sig- 
nificantly reduced due to the fact that the Mgll sam- 
ple size was increased 10-fold. Second, the number of 
DLAs in each bin has increased by more than a factor of 
3. Thus, the uncertainties in the low- and high-redshift 
data points are now comparable. Note that the statis- 
tics of the high-redshift data are also improved due to the 
inclusion of an SDSS DLA sample (Prochaska & Herbert- 
Fort 2004). Nevertheless, our basic conclusion from RT00 
has remained unchanged, namely, that the cosmologi- 
cal mass density of neutral gas remains roughly constant 
from z « 5 to z k 0.5. The drop in redshift number 
density from z = 5 to z = 2 along with a constant mass 
density in this range indicates that while the product 
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of galaxy cross section and comoving number density is 
declining, the mean column density per absorber is in- 
creasing. This is, again, consistent with the assembly of 
higher density clouds as galaxy formation proceeds. 

On the other hand, a constant cross section from z 1 
to z — along with a drop in mass density from z ~ 0.5 
to z — is indicative of star formation that depletes the 
highest column density gas while keeping the absorption 
cross section constant. This would in turn require that 
the column density distribution of DLAs change such 
that the ratio of high to low column densities decreases 
from low-redshift to z = 0. As we will see in the next 
section, the column density distribution does show some 
evidence for this. 

4.4. Column density distribution f(N) 

Figure 17 shows the normalized cumulative column 
density distribution (CDD) for the three redshift regimes. 
The dashed curve is the z = CDD from an analysis 
of an HI diameter-limited sample of local galaxies from 
Rao & Briggs (1993) while the dotted curve is the z — 
distribution derived by Ryan- Weber et al. (2003, hence- 
forth R-W03; 2005) from HIPASS data. The thick, solid 
curve is derived from the DLAs in Table 1 and the thin, 
solid curve is from the "total" sample of Prochaska & 
Herbert-Fort (2004). The change in the three CDDs 
with redshift is exactly what is expected based on the 
n DLA and Qdla results. Namely, that the low-redshift 
CDD shows a higher incidence of high column density 
systems than at high redshift presumably due to the as- 
sembly of gas as galaxy formation proceeds, followed by 
a decrease in the fraction of high column density sys- 
tems to z = 0, presumably due to the depletion of gas 
during star formation. Thus, at least qualitatively, the 
evolutionary behavior of Udla, &dla, and the CDD are 
entirely consistent with one another. A KS test shows 
that there is a 25% probability that the high- and low- 
redshift curves are drawn from the same population; this 
is significantly higher than what we observed in RT00, 
where the two samples had only a 2.8% probability of 
being drawn from the same population. However, the 
general trend that the low-redshift sample has a higher 
fraction of high column density system still remains. 

The absolute CDD can be determined using the equa- 
tion 

f(N,z) = n DLA (z) {1 + z)2 AN (6) 

where y(N, z) is the fraction of DLAs with column den- 
sities between N and N + AN at redshift z, and E(z) 
is as given in Equation 5. Figure 18 is a plot of the log 
of the absolute CDD function, log /(TV), as a function 
of log N (HI). The turnover with redshift is most ap- 
parent in the lowest and highest column density bins. 
We derive (3 = 1.4 ± 0.2 and f3 = 1.8 ± 0.1 at low 
and high redshift respectively, where the CDD is ex- 
pressed as f(N) = BN-P. At z = 0, R-W03 derive 
= 1.4 ± 0.2 for log N (HI) < 20.9 and (3 = 2.1 ± 0.9 
for log N(HI) > 20.9. The general form of the absolute 
CDD does not vary considerably with redshift, which in 
turn explains the roughly constant value of Qdla- The 
differences in the f(N) distributions are subtle, implying 
that the gas content in DLAs is not changing drastically. 
This is strong evidence that DLAs do not have high SFRs 
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Fig. 17. — The normalized cumulative column density distri- 
bution of DLAs for three redshift regimes. The top curve (red, 
thick solid line) includes the 41 low-redshift DLAs from Table 1 
at a median redshift of 0.95. The middle curve (blue, thin solid 
line) includes 163 high-redshift systems with mean redshift 2.94 
(Prochaska & Herbert- Fort 2004), and the bottom two curves are 
estimates at z = 0. The dashed curve is from Rao & Briggs (1993) 
and the dotted curve is from R-W03. 
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Fig. 18. — Absolute column density distribution (CDD) function 
for the three redshift regimes. The red, thick solid line is a least 
squares fit to the low redshift data points and has slope (3 = 1.4; 
the blue, thin solid line is a least squares fit to the high redshift 
data points with a slope of /3 = 1.8. The dotted line is the column 
density distribution derived by R-W03 from 21 cm HIPASS data 
of local galaxies. It has slope j3 = 1.4 for log N(H I) < 20.9 and 
P = 2.1 for log N (HI) > 20.9. The dashed line is from HI mea- 
surements of an optically-selected sample of local galaxies (Rao & 
Briggs 1993). The offset between the two z = curves arises from 
the different normalizations of the Hi-mass and optical luminosity 
functions, respectively. See text. 



and are, therefore, a different population of objects than 
those responsible for much of the observed luminosity in 
the high redshift universe (see also §5.1 and Hopkins et 
al. 2005). On the other hand, a non-evolving DLA popu- 
lation might be observed if the gas that is used up in star 
formation is replenished from the inter-galactic medium 
at a comparable rate. This possibility seems rather con- 
trived, and requires more proof than the current obser- 
vational evidence can provide. In §5, we discuss further 
evidence that suggests that DLAs and high SFR galax- 
ies (e.g., Lyman break galaxies) are mutually exclusive 
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populations. 

We'd like to draw the reader's attention to the z = 
curves in Figures 17 and 18. The general shape of the 
CDD at the present epoch is reproduced in both the 
RB93 and R-W03 estimates (Figure 17), with only a mi- 
nor difference in the column density where the slope of 
the distribution changes (Figure 18). This break in the 
CDD is a direct consequence of the exponential distribu- 
tion of N(HI) in disks, with the position of the break 
depending on the maximum column density in a face- 
on galaxy. The RB93 CDD was determined using HI 21 
cm maps of an optical diameter-limited sample of local 
galaxies normalized by the luminosity function of late- 
type galaxies locally (Rao 1994). On the other hand, the 
R-W03 estimate was derived from the HIPASS survey 
of HI in the local universe and the HI mass function of 
galaxies locally (Zwaan et al. 2003). The offset between 
the two curves (Figure 18) is a result of different normal- 
izations in the luminosity and HI mass functions of the 
two estimates, respectively. We now know that the local 
galaxy luminosity function (as determined prior to 1993) 
did not include gas-rich galaxies that occupy the low- 
luminosity tail of the galaxy luminosity function. The 
overall normalization might have been underestimated as 
well. This offset is also manifested in the values derived 
for n(z), since n(z) depends on the volume number den- 
sity of absorbers. In Figure 14, we used the most recent 
estimate of n(z = 0) = 0.045 ± 0.006 derived by Zwaan 
et al. (2005a). Ryan- Weber et al. (2005) derive a similar 
value; these are a factor of 3 higher than the RB93 esti- 
mate. It is also of interest to note that the local value of 
VIhi derived by RB93 is only ~ 30% smaller than more 
recent estimates. This means that while the bulk of the 
neutral gas at z = is in the more luminous galaxies, a 
significant fraction of the HI cross section is contributed 
by optically low-surface-brightness and dwarf galaxies. 
This conclusion is indeed borne out by the recent com- 
parison of HI in low- and high-surface brightness galaxies 
by Minchin et al. (2004) . In any case, it is now clear that 
the recent deep, large scale surveys of HI gas in the lo- 
cal universe have provided a better understanding of the 
distribution of HI at z = 0, allowing for more precise de- 
terminations of its statistical properties for comparison 
with quasar absorption line studies. 

5. DISCUSSION 

5.1. Selection Effects and Biases 

As with any survey, selection effects and biases need to 
be well understood in order to correctly interpret results. 
Here we raise some of the important ones that may affect 
our survey. 

1. We have determined DLA statistics under the as- 
sumption that all DLAs exhibit Mgll absorption, and 
therefore, that DLAs form a subset of Mgll absorbers. At 
the redshifts probed by our UV surveys, 0.11 < z < 1.65, 
we find that there is little chance of encountering a DLA 
unless Wo A2796 > 0.6 A. Since our Mgll sample includes 
systems with W A2796 > 0.3 A, we believe this result to 
be fairly secure. An exception may occur in the rare 
case where the DLA sightline passes through a single 
cloud. Its velocity width, i.e., b parameter and there- 
fore, W^ 2796 , would then be small, perhaps even smaller 
than in the DLA towards 3C 286 (see §4.1.1). The DLA 



towards the D component of the Cloverleaf gravitation- 
ally lensed quasar might be an example of such a case 
(Monier et al. 2005). The N(HI) = 2 x 10 20 cm" 2 , 
z = 1.49 DLA is not detected in the three brighter com- 
ponents of this quadruply lensed quasar. A composite 
Keck spectrum of all four components combined does 
not detect Mgll down to a 3er limit of Wg A2796 = 0.06 A. 
Given the relative brightnesses of the four components, 
any Mgll absorption towards component D could be di- 
luted by a factor of 5 - 10 in the composite spectrum. 
Thus, the absence of metal lines in this DLA need not 
be an indication of unpolluted gas, but instead, of low 
velocity dispersion gas that might only be detected with 
high resolution spectra of component D. 

Although the highest redshift DLAs have not been 
shown to have Mgll absorption because the Mgll doublet 
is shifted into the infrared, metals have been detected 
in DLAs at redshifts as high as z — 3.9 (Prochaska et 
al. 2003a) and are, therefore, expected to also include 
Mgll. Apart from testing the DLA-Mgll connection at 
high redshift and exploring any evolution, assembling a 
near-infrared Mgll sample with follow-up optical spec- 
troscopy to search for DLAs would be important for com- 
parison with blind optical DLA surveys. In addition, 
Nestor et al. (2003) and Turnshck et' al. (2005) have 
shown that a positive correlation exists between Wq 2796 
and neutral gas metallicity when an ensemble average of 
strong (W A2796 > 0.6 A) SDSS Mgll absorbers is consid- 
ered. Since the REW of saturated Mgll lines is an indi- 
cation of gas velocity spread, this correlation indicates 
a metallicity-kinematics relation for the average Mgll 
absorber. The evolution of this relationship to higher 
redshift will also provide important constraints on CDM 
simulations of galaxy and structure formation. 

Ultimately, our study is based on the premise that 
all DLAs have Mgll absorption, and unless significant 
counter examples are found, this assumption is now on 
fairly firm ground. 

2. Figure 4 shows that the fraction of DLAs in a Mgll 
sample is a function of Wq 2796 , rising from a fraction near 
16% just above the threshold value of W A2796 = 0.6 A to 
about 65% at the highest values near W^ 2796 — 3 A. At 
present these should be considered approximate fractions 
since the presence of FeIIA2600 has been used to increase 
the probability of finding a DLA. However, for reasons 
that are not yet clear, the Fell criterion does not affect 
the mean HI column density as a function of Wq 2796 (see 
Figure 7), and so the Fell inclusion effect is probably 
not significant for our sample. Nevertheless, the Wq 2796 
dependence of the DLA fraction will introduce a bias in 
nDLA unless the observed sample's Wq 2796 distribution 
matches the true distribution. We have accounted for 
this by carefully defining our samples as described in §2, 
and by making use of Equations 1,2, and 3 to calculate 

nDLA- 

3. The degree to which N(HI) may be biased by 
Wq 2796 can be seen by examining Figure 5. The mean HI 
column density of identified DLAs is N(HI) « 2 x 10 21 
cm" 2 when 0.6 A < VF A2796 < 1.2 A, but it seems to de- 
crease by a factor of « 4 at W^ 2796 « 3 A. However, 
inspection of Figure 2 suggests that this trend is not 
particularly tight nor is it well established for DLAs by 
themselves. On the other hand, if one considers all the 
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points in Figure 2, it is clear that the N(HI) distribu- 
tion changes for different Wq 2796 intervals. It is inter- 
esting that in Mgll-sclected surveys for DLAs, the de- 
termination of the cosmological mass density of neutral 
gas, QdlAi has (so far) not revealed any dependency 
on Wq 2796 selection. This is because in our sample, 
the increased probability (by a factor of « 4) of find- 
ing a DLA at the largest Wq 2796 values is approximately 
compensated for by the corresponding decrease in mean 
HI column density (by a factor of w 4) at the largest 
Wq 2796 values. It is worth pointing out that although 
the Mgll selection criteria lead to reasonably complete 
samples of DLAs, incompleteness must set in at HI col- 
umn densities in the sub-DLA regime because systems 
with Wq 2796 < 0.3 A can have sub-DLA HI column densi- 
ties. Therefore, only the N(HI) distribution in the DLA 
regime can be reliably considered with the available data. 

4. Hopkins ct al. (2005; see also Rao 2005) have 
discussed the question of whether the observed popu- 
lation of DLAs can account for the observed SFH of 
the universe from low to high redshift. By applying the 
Kennicutt (1998) formulation of the Schmidt law to the 
properties of the currently observed population of DLAs 
they find that the DLAs cannot account for the cosmic 
SFR density inferred from the luminosity density of high- 
rcdshift galaxies (see figure 2 in Hopkins et al. 2005). An 
even larger discrepancy occurs when one compares DLA 
metallicities to the metallicities expected on the basis of 
the cosmic SFR (see figure 4 in Hopkins et al. 2005 and 
figure 13 in Rao 2005). One way to avoid this discrepancy 
is to postulate that the Mgll and blind DLA surveys are 
not yet large enough to include absorbers with very small 
individual cross sections that nevertheless may dominate 
the cosmic SFR and be the main reservoirs for the met- 
als as well. Indeed, these star forming regions will be 
rich in molecular gas, the direct fuel for star formation, 
but with HI column densities that may exceed the ob- 
served DLA regime. Kennicutt (1998) points out that 
in normal disks star formation generally takes place in 
regions that contain 1 - 100 M Q pc -2 (i.e., « 10 20 - 10 22 
atoms or molecules cm -2 ), whereas the more rare (and 
smaller) star burst regions contain 10 2 — 10 5 M pc~ 2 
(i.e., w 10 22 — 10 25 atoms or molecules cm" 2 ). For ex- 
ample, an absorber with a size of about 100 pc, com- 
parable to giant molecular clouds (GMCs), has a cross 
section that is w 10 4 times smaller than known DLAs, 
which typically have effective radii of « 10 kpc (Monicr 
et al. 2005). Assuming that there are on the order of 
10 GMCs per galaxy, the total cross section per unit 
volume, i.e., interception probability, for these very high 
column density gas systems would be on the order of 
10 3 times smaller. This means that 10 3 DLAs need to 
be detected in order to find one very high column den- 
sity system. With the SDSS, we are getting close, but 
are not quite there yet. A one in a thousand system 
with N(HI + H 2 ) = 10 24 cm" 2 would increase the SFR 
density of DLAs by more than a factor of 2, and bring 
the DLA SFR density into agreement with the luminous 
SFR density. Searches for molecular gas in DLAs have 
resulted in only a handful of detections. Moreover, the 
molecular gas fraction in the few DLAs with H2 detec- 
tions is very small (e.g. Ledoux et al. 2003), and is 
consistent with the idea that the known sample of DLAs 



does not trace the majority of the star forming gas in the 
universe. It therefore seems reasonable to conclude that 
most of the neutral and molecular gas mass has so far 
been missed in DLA surveys. 

However, the possibility that these very high column 
density gas systems are being missed by DLA surveys 
may not only be due to their small gas cross sections, but 
also because they are likely to be very dusty. Ledoux et 
al. (2003) find that the DLAs in which H2 is detected 
have among the highest metallicities and the highest de- 
pletion factors, hinting at the possibility of much higher 
depletions in much higher column density molecular gas 
clouds. Although radio loud quasar surveys for DLAs 
have not revealed any significant dust bias in optical sur- 
veys (Ellison et al. 2001; 2004; Akerman et al. 2005), 
the radio loud quasar surveys for DLAs may themselves 
suffer from the small cross section selection effect. Not 
enough radio loud quasars have yet been surveyed to find 
the putative one in a thousand very high column density 
system. But if significant mass has been missed due to 
small total cross section for star forming regions, whether 
or not these high-gas-mass regions will be found once 
sample sizes are much larger is unclear. Substantial dust- 
induced reddening may prevent complete samples from 
ever being discovered via optical quasar absorption-line 
spectroscopy. 

In this regard it is interesting that Gardner et al. 
(1997) found in their CDM simulations that depletion of 
the gas supply by star formation only affected the DLA 
statistics at z > 2 for N(HI) > 10 22 atoms cm" 2 (i.e., in 
a regime where DLAs have not been found), even though 
roughly half of the cold collapsed gas was converted to 
stars by z = 2. 

5. Gravitational lensing has the opposite effect on 
DLA surveys. Magnification by DLA galaxies could 
brighten background quasars, and preferentially include 
them in magnitude-limited samples. Le Brun et al. 
(2000), with HST imaging observations, showed no ev- 
idence for multiple images of background quasars and 
concluded that the quasars were magnified by at most 
0.3 magnitudes. In addition, Ellison et al. (2004) and 
Peroux et al. (2004) using statistical tests on low redshift 
Mgll and DLA samples, showed that lensing bias is a mi- 
nor effect. More recently, using the SDSS Mgll survey 
results of Nestor (2004), Menard et al. (2005) show that 
quasars behind strong Mgll absorbers, of which DLAs 
are a subset, show little magnification bias, and that 
its effect on Qdla at low redshift is negligible (see also 
Menard 2005) . It is also unlikely that the lowest redshift 
points that we derived from our HST-UV data (Figure 
16) are affected by lensing bias. This is because the DLAs 
with the highest HI column densities at z « 0.5 arise in 
dwarf galaxies (Rao et al. 2003), and consequently, do 
not have the mass required to produce significant mag- 
nification. 

5.2. Interpretation of the Statistical Results on DLAs 

As discussed in §4, the evolutionary behavior oin^LA, 
Qdla, arid the CDD are, at least qualitatively, consis- 
tent with one another in terms of a simple galaxy for- 
mation scenario. We find evidence for a rapid decline of 
n DLA (by a factor of 2) from z = 5 to z « 2 followed by 
no evolution down to z = 0. For comparison, the evo- 
lution of Lya forest lines with log N (HI) > 14, which 
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Fig. 19. — Redshift number density of DLAs as a function of 
redshift. Symbols are the same as in Figure 14 with the addition 
of open triangles, which are derived from the Wq 27 ® & > 0.6 A Mgll 
redshift number density and assuming that the fraction of DLAs 
in these Mgll systems is constant at 22%. The errors are there- 
fore indicative of statistical errors in the Mgll sample alone. The 
open star at z = 0.06 is similarly derived from the HST Mgll sam- 
ple of Churchill (2001). Including errors in the DLA fraction will 
systematically move the data up or down by ~ 25%. 
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Fig. 20. — Cosmological mass density of neutral gas in DLAs 
as a function of redshift. Symbols are the same as in Figure 16 
with the addition of open triangles, which are derived from the 
IVq 12796 > 0.6 A Mgll redshift number density, and assuming that 
the fraction of DLAs in these Mgll systems is constant at 22% 
and that their HI column density is constant at 1.16 X 10 21 cm -2 . 
The errors are therefore indicative of statistical errors in the Mgll 
sample alone. The open star at z = 0.06 is similarly derived from 
the HST Mgll sample of Churchill (2001). Including errors in the 
DLA fraction and N(HI) will systematically move the data up or 
down by ~ 0.1 dex. We see, for the first time, possible evidence of 
a decline in £Idla from z ?s 0.5 to z = 0. 



have been shown to be associated with the same large 
scale structure that traces galaxies (Tripp, Lu, & Savage 
1998; Dave et al. 1999; Penton, Stocke, & Shull 2002), 
also slows down dramatically near z«1.5 (e.g., Kim et 
al. 2002, Weymann et al. 1998). Thus, both the DLAs 
and the higher N(HI) Lya forest appear to follow simi- 
lar evolutionary histories consistent with the collapse and 
assembly of baryonic structures near z ~ 1.5 or 2. The 
near constant value of floLA during the phase of rapid 
evolution in uola implies an increase in the HI column 
densities of individual clouds. The observed evolution 
in the CDD of the DLAs, although mild, is evidence for 
this. The subsequent drop in Qdla down to z = along 
with an unevolving ji£> la is indicative of star formation 
that depletes gas while keeping the absorption cross sec- 
tion constant. The change in slope of the CDD from low 
redshift to z = 0, i.e., the decrease in the ratio of high to 
low column densities, is again consistent with this sce- 
nario. Further details on the evolution of HI from low 
redshift to the present epoch can be studied only when 
the sample of low-redshift DLAs becomes large enough to 
split the < z < 1.65 redshift interval, without compro- 
mising on the uncertainties, into finer than the current 
two bins, and now we consider this possibility by adopt- 
ing some reasonable assumptions. We have shown 
that a survey of Mgll systems with Wq 2796 > 0.6 A is a 
reliable tracer of DLAs, and can be used to determine 
DLA statistics. The two hola data points at low red- 
shift shown in Figure 14 are 20 (± 5)% and 24 (± 6)% of 
the corresponding Mgll redshift number density values 
for Wo A2796 > 0.6 A derived by Nestor (2004) and NTR05, 
at z = 0.6 and z = 1.2 respectively. By assuming that 
the DLA fraction in a Wq 2796 > 0.6 A sample is constant 
over the entire redshift interval 0.1 < z < 1.65, tidla 
can be estimated in much smaller redshift bins from the 



statistics of any Mgll survey sample without a UV sur- 
vey for DLAs. The systematic uncertainty in udla will 
then be primarily limited only by the precision to which 
the DLA fraction is known. Similarly, the systematic 
uncertainty in Qdla will primarily be limited by how 
accurately the DLA HI column density is known. For 
Wq 2796 > 0.6 A, we found the mean DLA column den- 
sity to be (N{HI)) = (1.16 ±0.20) x 10 21 cm" 2 . Recall 
that all the DLAs in our sample have VF A2796 > 0.6 A. As- 
suming a constant DLA fraction of 22% and a constant 
DLA HI column density of 1.16 x 10 21 cm -2 , we can es- 
timate n DLA and Q, DLA from the Mgll VF A2796 > 0.6 A 
redshift number density as follows (see Equations 1 and 
4): 

n D LA{z) = 0.22 x n Mg ii{z) (7) 

and 



n nLA (z) = Mmgg ° 0-22 xn Mg ii(z) x 1.16 x 10 21 



E{z) 



CPc 



(1 + z) 2 ' 
(8) 

These data points are shown in Figures 19 and 20 as open 
triangles. Only the error in riMgii is propagated through 
to show the statistical uncertainty in these data. The er- 
rors associated with the DLA fraction and DLA column 
density are systematic and will affect all of these data 
points equally, moving them uniformly up or down by 
~ 25% in the case of udla and ~ 0.1 dex for \og£lDLA- 
We also show tidla and Qdla inferred using the red- 
shift number density derived by Churchill (2001) for Mgll 
systems detected in HST spectra, riMgii(z = 0.06) = 
0.221^09! again, only the error in riMgii has been prop- 
agated. We now see for the first time that, under the as- 
sumption of constant DLA fraction and HI column den- 
sity suggested by our current Mgll-DLA sample, there 
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Fig. 21. — Redshift number density of DLAs as a function of 
cosmic time with to being the current epoch. Symbols are the 
same as in Figure 19. 
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Fig. 22. — Cosmological mass density of neutral gas in DLAs as 
a function of cosmic time with to being the current epoch. Symbols 
are the same as in Figure 20. 



may be evidence of a decreasing trend in Hdla from 
z = 0.5 to z = 0. There also appears to be a dip in 
Qdla near z = 2, albeit within la, that shows up in 
both the high- and low-redshift data. It should be noted 
that the highest of the low redshift data points comes 
from the red end of SDSS spectra, and suffers from low 
signal-to-noise ratio due to the presence of atmospheric 
absorption. Similarly, the lowest of the high redshift data 
points comes from the blue end of SDSS spectra, and also 
suffers from low signal-to-noise ratio. Whether this la 
effect will persist with better quality data remains to be 
seen; but if real, will be a challenge for galaxy forma- 
tion models to explain. These two figures illustrate 
our method for determining DLA statistics. Using our 
sample of 197 Mgll systems with follow-up UV spectra, 
we have demonstrated that, to first order, the DLA frac- 
tion in a W A2796 > 0.6 A Mgll sample and the DLA 
HI column density are constant. Using this assumption, 
we have shown that details in the evolution of Mgll sys- 
tems can reveal details in the neutral gas evolution. The 
data point at z = 0.28 was derived from the MMT sur- 
vey for low redshift Mgll systems (Nestor 2004; Nestor, 
Turnshek, & Rao, in preparation), and is in a redshift 
regime inaccessible by the SDSS. An even larger survey 
for Mgll systems at 0.11 < z < 0.36 is clearly needed in 
order to understand the evolution of Qdla hi a redshift 
regime where most of the evolution appears to be taking 
place. Similarly, a high-redshift, near infrared survey for 
Mgll systems could extend this method into the optical 



regime, and any evolution in the Mgll-DLA relationship 
could be studied. In order to underscore the importance 
of pursuing this work in the future, we show uola and 
Qdla as functions of cosmic time in Figures 21 and 22, 
respectively. 

6. SUMMARY 

We have presented statistical results on UV surveys 
for low-redshift (z < 1.65) DLAs with N(HI) > 2 x 10 20 
cm~ 2 using the largest sample of UV-detected DLAs ever 
assembled. The DLAs were found by targeting QSOs 
with Mgll systems identified optically in the redshift 
range 0.11 < z < 1.65. In total, UV observations of the 
Lya absorption line in 197 Mgll systems with Wq 2796 
> 0.3 A have been obtained. This is an efficient and ef- 
fective way to find DLAs because, in the absence of Mgll, 
the system evidently has no chance of being a DLA. This 
sample contains 41 DLAs, all of which have Wq 2796 > 0.6 
A. Our main findings can be summarized as follows: 

1. To a high level of completeness, DLAs can be stud- 
ied through follow-up observations of strong Mgll ab- 
sorbers. In particular, Figure 4 shows that for our sample 
the probability of a Mgll system being a DLA is P « 
for TU A2796 < 0.6 A and P w 0.16 + 0.18(Wo A2796 -0.6) for 
0.6 < Wo A2796 < 3.3 A. A Mgll absorber must generally 
have 1 < W A2796 /W A2600 < 2 and 1U A2852 > O.lA to be a 
DLA (see §4.1). 

2. UV spectroscopy, almost exclusively with HST, en- 
abled us to measure or place limits on N(HI) for each 
of the 197 systems studied. For Mgll systems with 0.3 
A < 1U A2796 < 0.6 A, (N(HI)) = (9.7 ± 2.2) x 10 18 
cm" 2 , while for systems with W A2796 > 0.6 A, (N(HI)) = 
(3.5 ± 0.7) x 10 20 cm~ 2 . This is basically a step func- 
tion (see Figure 4) , with a factor of ss 36 change in mean 
HI column density near W A2796 « 0.6 A. Since the Mgll 
absorption lines are saturated at Wq 2796 > 0.6 A, there 
is evidently a threshold in kinematic velocity spread be- 
low which it is highly unlikely to encounter high column 
density neutral DLA gas. 

3. Above IU A2796 = 0.6 A, the mean HI column density 
of a sample of Mgll absorbers is found to be constant 
with increasing Wq 2796 . However, owing to the increase 
in probability of finding a DLA with increasing Wq 2796 , 
the mean HI column density of Mgll absorbers that are 
DLAs is found to decrease by about a factor of four with 
increasing W A2796 , from W^ 2796 ^ 0.6 A to W A2796 ss 3.5 
A. Improved statistics are needed to study this effect 
owing to the large scatter in the W / A2796 versus N(HI) 
plane for W A2796 > 0.6 A. 

4. By combining results at all redshifts, including 
21 cm emission surveys at z = 0, we find that the 
DLA incidence per unit redshift can be parameterized 
as n DLA (z) = n (l + z) 7 where n = 0.044 ± 0.005 
and 7 = 1.27 ± 0.11. In the standard "737" cosmology 
this indicates no evolution in the product of neutral gas 
cross section times comoving number density at redshifts 
z < 2, but from z ss 5 to z « 2 there is a decrease of a 
factor of ss 2 in this quantity relative to the no evolution 
prediction (Figure 19). This decline happens relatively 
rapidly, in a time span that corresponds to < 1.5 Gyrs. 

5. The cosmological mass density of neutral gas due to 
DLAs, QdlAi follows a completely different evolutionary 
pattern. It remains relatively constant in the redshift in- 
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terval 0.5 < z < 5.0, with FIdla ~ 10~ 3 , but then it 
declines by a factor of « 2 between z « 0.5 and z — 
(Figure 20). This drop in neutral gas takes place during 
the last w 5 Gyrs of the history of the Universe. However, 
due to possible selection effects which are biased against 
finding regions with very high column densities because 
the product of their gas cross section and comoving num- 
ber density is small, it is important to realize that the 
neutral gas component as traced by the DLAs may not 
include all of the neutral and molecular gas involved in 
star formation (Hopkins et al. 2005; Rao 2005; Turnshck 
et al. 2005). 

6. Consistent with the uy)la{z) and £Idla(z) results, 
the HI CDD at (z) « 1 shows a higher incidence of high 
column density systems than at (z) w 3. This presum- 
ably represents a build up of neutral mass concentrations. 
By z — the higher incidence of high N(HI) systems 
seen at (z) w 1 has disappeared, presumably due to the 
depletion of gas during star formation. 

7. In the absence of future QSO absorption-line sur- 
veys that aim to identify DLAs and measure their N(HI) 
in UV spectra, more detailed studies that lead to a bet- 
ter understanding of the strong Mgll systems may hold 
promise for reaching a better determination of the prop- 



erties of the neutral gas phase of the universe at z < 1.65 
(e.g., Figures 19 - 22). 
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n Q/i 9K 

U.o4Z0 
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± 


0.038 
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1.361 


± 
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1 
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1 
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± 
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± 
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4 
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Zl.oz_ Q7 
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± 
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± 
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4 
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± 


0.047 
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4 
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1 
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5 
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1 
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± 


0.05 


0.68 


± 


0.05 
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Z 
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9 

z 


HQ 

uy 
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y 


1 
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0.862 


± 


0.074 
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± 


0.070 
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u.uyo it u.uoo 


A 
4 
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1 
1 
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18 


2 


1 
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0.6322 


1.140 ± 


101 


1.579 


± 


0.087 


1.216 


± 


0.093 


0.568 ± 0.091 


4 


19 95 + nnl 


4 


11 












0.7087 


0.888 ± 


091 


1.210 


± 


0.066 


1.008 


± 


0.068 


0.718 ± 0.088 


4 


2o.o4i8;8I 


4 


11 


1032+0003 


18 


9 


1 
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1.0168 


0.697 ± 


119 


1.919 


± 


0.128 


1.650 


± 


0.130 


-0.083 ± 0.160 


4 


i9.ool8;? 


4 


11 


1035-276 


19 





2 


168 


0.8242 


0.52 ± 


02 


1.08 


± 


0.02 


0.87 


± 


0.02 


<0.6 


26 


i8.8it8:IS 


1 


06 


1037+0028 


18 


4 


1 


733 


1.4244 


1.767 ± 


081 


2.563 


± 


0.031 


2.171 


± 


0.031 


0.303 ± 0.042 


4 


20.04181° 


4 


11 


1038+064 


16 


7 


1 


265 


0.4416 


<0.2 




0.66 


± 


0.05 


0.57 


± 


0.04 


<0.2 


2 


18-3018:38 


1 


AR 


1040+123 


17 


3 


1 
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0.6591 


<0.2 




0.58 


± 


0.10 


0.42 


± 


0.10 


<0.2 


14 
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1 


AR 


1047-0047 


18 
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0.5727 


0.765 ± 
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1.063 


± 


0.117 


0.697 


± 


0.094 


0.055 ± 0.096 


4 


19-3618;^ 


4 


11 


1048+0032 


18 


6 


1 


649 


0.7203 


1.252 ± 


077 


1.878 


± 


0.063 


1.636 


± 


0.070 


0.417 ± 0.085 


4 


18.78l8il 


4 
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TABLE 1 — Continued 



QSO 




z a.bs 


^2600 (A) 


W X2796 (A) 


W X2803 (A) 


w \2«a { A) 


Refs 


log N(HI) b Sel c 


UV sp d 






Mgll 


Fell 


Mgll 


Mgll 


Mgl 




Crit 


source 



1049+616 
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.421 


0.2255 


<0.4 




0.51 ± 0.03 


0.56 ± 0.03 


<0.1 


18 
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1 


06 
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0.5604 


2.132 ± 0. 


117 


2.273 ± 0.089 


2.623 ± 0.083 
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AR 
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AR 


1254+047 
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1 


.018 
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0.40 ± 0. 


05 


0.46 ± 0.04 


0.38 ± 0.04 


<0.2 
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1 


AR 
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1 
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1 


AR 
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1 
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0.7160 
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4 


20.54 


4 


11 
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17.5 


1 
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0.57 


0.53 
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3 
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1 
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0.16 


3 
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1 


06 
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1 
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32,33 
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1 


AR 
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1 


937 
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<0.3 
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2 
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1 


06 
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1 
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<0.1 
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0.17 ± 0.04 


<0.3 


2 


19.081°°! 


1 


AR 
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18.8 


1 


.714 
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a For systems other than those from Nestor (2004), upper limits are either given by the authors or estimated by us from the published spectrum. If the line is part of a blend, then 
the rest equivalent width entered is an upper limit. The upper limit is a la upper limit if estimated by us from the published spectrum. In some cases, the absorption line (A2600 or 
A2852) was clearly visible in the spectrum but was not identified by the authors since the rest equivalent width of the line did not meet their detection criterion, usually 5<r. In these 
cases, we estimated the equivalent width of the line from the spectrum and tabulated the measurement as an upper limit. The trends that we have established are evident from the 
measurements of rest equivalent widths by the authors of the published spectra; our upper limit measurements only serve to provide a certain degree of completeness to our sample 
and, in fact, do not indicate or contribute to any trends by themselves. We note that the quasar 0151+045 was erroneously included in Table 4 of RT00. The Mgll system towards this 
quasar was detected after the galaxy-quasar pair was known and is, therefore, a biased system. In addition, the z a \y S = 0.213 system towards 1148+386 and the z a \y S = 0.1634 system 
towards 1704+608 were flagged as doubtful systems by Boisse et al. (1992). Finally, the nature of the z a i s = 1.3284 system towards 1331+170 was deemed inconclusive from the low 
signal-to-noise ratio IUE spectrum. Therefore, these four were eliminated from our current Mgll sample. 

b In two cases (the z a t, s = 1.6101 system towards 1329+412 and the z a i, s = 1.2528 system towards 1821+107), the HI column density could not be determined because the Lya and 
higher order Lyman lines could not be fit with a unique value of N(HI) under the assumption that the lines are damped. A large b value (greater than 100 km s _1 ) was necessary in 
both cases. Thus, it is clear that these systems are not DLAs. These are flagged as "high b" systems. 

c Flag for Wq 2796 selection criterion used in the determination of «_d_la- Sec §4. 

d Source for UV spectrum from which Lya information was obtained. AR: HST or IUE archive, 06: HST-cycle 6 program 6577, P.I. Rao, 09: HST-cycle 9 program 8569, P.I. Rao, 
11: FST-cycle 11 program 9382, P.I. Rao 



